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ABSTRACT 

We use a large set of state-of-the-art cosmological N-body simulations [512 3 particles] 
to study the intrinsic ellipticity correlation functions of halos. With the simulations 
of different resolutions, we find that the ellipticity correlations converge once the ha- 
los have more than 160 members. For halos with fewer members, the correlations are 
underestimated, and the underestimation amounts to a factor of 2 when the halos 
have only 20 particles. After correcting for the resolution effects, we show that the 
ellipticity correlations of halos in the bigger box (L — 300 /i _1 Mpc) agree very well 
with those obtained in the smaller box (L = 100 /i _1 Mpc). Combining these results 
from the different simulation boxes, we present accurate fitting formulae for the ellip- 
ticity correlation function cn(r) and for the projected correlation functions Sn(r p ) 
and ^22{fp) over three orders of magnitude in halo mass. The latter two functions 
are useful for predicting the contribution of the intrinsic correlations to deep lensing 
surveys. With reasonable assumptions for the redshift distribution of galaxies and for 
the mass of galaxies, we find that the intrinsic ellipticity correlation can contribute 
significantly not only to shallow surveys but also to deep surveys. Our results indicate 
that previous similar studies significantly underestimated this contribution for their 
limited simulation resolutions. 

Key words: Cosmology: theory — cosmology: observations — gravitational lensing 
— large-scale structure of universe — dark matter 



1 INTRODUCTION 

Inhomogeneities of matter distribution in the Universe dis- 
tort the images of distant galaxies gravitationally, a phe- 
nomenon called gravitational lensing. The lensing effect in- 
duces an ellipticity-ellipticity correlation of the galaxies on 
large scales, which is observable and can be used as a power- 
ful tool to probe the large-scale dark matter distribution in 
the Universe (Miralda-Escude 1991; Blandford et al. 1991; 
Kaiser 1992; see Bartelmann & Schneider 2001 and Mellier 
1999 for reviews). Several groups have already detected the 
ellipticity correlation on scales from 0.4 to 30 arc-minutes 
for faint galaxies (Bacon et al. 2000, 2002; Hoekstra et al. 
2002a,b; Kaiser et al. 2000; Maoli et al. 2001,; Rhodes et 
al. 2000, 2001; van Waerbeke et al. 2000, 2001; Wittman et 
al. 2000) . If the source galaxies are randomly oriented with- 
out the lensing effect, that is, the intrinsic ellipticity corre- 
lation of the galaxies is negligible, the observed ellipticity 
correlation implies that the parameter (3 = Qq 6 o-$ is about 
0.6, where f2o is the density parameter, as is the current rms 
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linear density fluctuation in a top-hat sphere of 8 ft 1 Mpc, 
and h is the Hubble constant in units of 100km s _1 Mpc _1 . 

There are however evidences both from theory and from 
observations that the shapes of galaxies are intrinsically cor- 
related. In the theory, the large scale tidal field is expected to 
induce large-scale correlations in galaxy spins and in galaxy 
shapes (Lee & Pen 2000, 2001; Croft & Metzler 2000, here- 
after CM00; Heavens et al. 2000, hereafter HRH00; Catelan 
et al. 2001; Catelan & Porciani 2001; Hui & Zhang 2002; 
Porciani et al. 2002). It is recently claimed that a large-scale 
alignment of galaxy spins has been detected in nearby galaxy 
catalogs with a high confidence (Pen et al. 2000, Lee & Pen 
2002; Brown et al. 2002; Plionis 1994 for cluster shapes). 
While the intrinsic ellipticity correlation may be separated 
from the weak lensing signal in observations through mea- 
suring the .E-mode and the B-mode correlations of the el- 
lipticity (Crittenden et al. 2001, Mackey et al. 2002), apply- 
ing this technique needs an accurate relation between the 
_E-mode and B-mode correlations (Crittenden et al. 2001, 
2002; Pen et al. 2002; Schneider et al. 2002; Hoekstra et al. 
2002b). In this aspect, the current situations are far from 
satisfactory, since there are still considerable uncertainties 
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both in the theoretical predictions and in the observational 
measurements of the intrinsic ellipticity correlations. 

Using the numerical simulations of 256 3 particles re- 
leased by the Virgo Consortium, HRHOO and CMOO mea- 
sured the ellipticity correlation for dark matter halos and 
discussed their results in the context of weak lensing mea- 
surement^ Assuming that the galaxies have the same in- 
trinsic correlations as their host halos, they found that the 
intrinsic correlation of galaxy ellipticity could dominate over 
the lensing signal in shallow lensing surveys, and could con- 
tribute a non- negligible signal, as high as 20 percent, to 
deep lensing surveys like those that reported detections of 
the weak lensing effects. Their findings already have pro- 
found implications for interpreting the weak lensing exper- 
iments, but quantitatively speaking, there are still signifi- 
cant uncertainties in their results. For examples, CMOO com- 
pared the ellipticity correlation cn(r) for halos of a mass 
> 1.4 x 10 12 /i _1 Mq between two simulations of boxsizes 
L = 141/i -1 Mpc and 240 /i -1 Mpc, and found that cii(r) 
in the higher resolution simulation (L — 141 ft _1 Mpc) is 
2 ~ 3 times higher. The ellipticity correlation cn(r) found 
in CMOO is also a factor of 2 or more higher than that found 
in HRHOO. Therefore, it is unclear how their selections of 
halos and simulation resolutions have affected their results. 

In this Letter, we will first quantify how the simulation 
resolution affects the determination of the ellipticity corre- 
lation using a set of state-of-the-art cosmological N-body 
simulations of 512 3 particles (Jing & Suto 2002) and a set 
of lower-resolution simulations of 256 3 particles (Jing 1998; 
Jing & Suto 1998). We will show that the minimum particle 
number iV m i n for resolving the halo shapes in simulations is 
160, and the ellipticity correlation functions of halos with 
more than iV m i n particles are well converged. We will also 
show that the ellipticity correlation is underestimated by 
a factor 2 when halos of about twenty particles are used. 
Therefore, the predictions of CMOO and HRHOO for weak 
lensing surveys are significantly underestimated, since the 
halos of ten (HRHOO) or twenty (CMOO) particles were in- 
cluded in their studies (in order to have a sufficient number 
of halos for their analyses) . Our results imply that the intrin- 
sic ellipticity correlation may have a significant contribution 
even to deep weak lensing surveys. Our simulations are also 
large enough for accurately measuring the scale- and mass- 
dependences of the ellipticity correlation functions. Based 
on the simulation data, accurate fitting formulae are pre- 
sented for these functions. These formulae can be used to 
predict the intrinsic ellipticity signals, including i5-mode and 
B-mode contributions, in large-scale lensing surveys (Crit- 
tenden et al. 2001, 2002; Schneider et al. 2001; Pen et al. 
2002). 

In the next section, we will measure the ellipticity cor- 
relation functions in a set of N-body simulations, and will 
do the convergence test. In Section 3, we will present the 
projected ellipticity correlation functions, which are useful 
for predicting the intrinsic ellipticity contribution in weak 
lensing surveys once the radial (or redshift) distribution of 
source galaxies is known. Our results are summarized and 
discussed in Section 4. 



HRHOO also examined the correlation of the halo spins. 



Table 1. Cosmological N-body simulations 



Llh-i-Mpc] 


N p 




m p [/i- 1 M ] 


realizations 


100 


512 3 


0.9 


6.2 x 10 8 


4 


300 


512 3 


0.9 


1.7 x 10 10 


4 


100 


256 3 


1.0 


5.0 x 10 9 


3 



2 THE ELLIPTICITY CORRELATION 
FUNCTIONS OF HALOS 

We use a set of cosmological N-body simulations of 512 3 par- 
ticles that are listed in Table 1. The cosmological model is 
a currently popular flat low-density model with the density 
parameter flo = 0.3 and the cosmological constant Ao = 0.7 
(LCDM). The shape parameter T = Qoh and the ampli- 
tude erg of the linear density power spectrum respectively 
are 0.2 and 0.9. Two different boxsizes L = 100/i _1 Mpc 
and L — 300fe _1 Mpc are adopted, and the particle mass 
m p is 6.2 x 10 8 /i~ 1 M Q and 1.7 x lO lo /i -1 M respectively 
in these two cases. Thus, halos of galactic sizes are well re- 
solved in both types of simulations. The simulation data 
were generated on the VPP5000 Fujitsu supercomputer of 
the National Astronomical Observatory of Japan with our 
vectorized-parallel P 3 M code (Jing & Suto 2002), and those 
of L — 100 /i^ 1 Mpc have already been used by Jing & Suto 
(2002) to derive a tri-axial model for density profiles of halos 
which is significantly more accurate than the conventional 
spherical model. There are 4 realizations for each boxsize. 
For comparison, we use three realizations of our previous N- 
body simulations for the LCDM model that contains 256 3 
particles (Jing & Suto 1998; Jing 1998). These simulations 
have the same initial fluctuation phases as the first three 
realizations of 512 3 particles and L = 100/i _1 Mpc, but the 
normalization ag, = 1-0 in the lower-resolution simulations 
is slightly higher. We refer readers to Jing & Suto (2002) for 
complementary information about the simulations. 

Dark matter halos are identified with the friends-of- 
friends method (FOF). A linking length b equal to 0.1 of 
the mean particle separation is adopted. This linking length 
is smaller than the nominal value 0.2. The regions selected 
b — 0.2 have a mean density contrast about 180 (Porciani 
et al. 2002), which approximately correspond to virialized 
halos. The objects identified with b = 0.1 have a mean over- 
density 8 times higher, and they are actually the inner re- 
gions of virialized halos (i.e. regions within about one third 
of the virial radius). Since we are interested in the elliptic- 
ity correlation of galaxies and the galaxies lie in the inner 
regions of halos, b — 0.1 might be more appropriate for our 
study (see also CMOO). The simulation outputs at redshift 
z = 1 are used, for galaxies in deep lensing surveys are at 
this redshift. We measure the ellipticity correlations that are 
defined (e.g., Miralda-Escude 1991; HRHOO; CMOO) as 

ey(r) =< ei (x)e f (x + r) > (1) 

where i and j runs over 1 and 2, and r is the three dimen- 
sional vector connecting a pair of halos. The ellipticity ei 
and ei for the halos are defined in a projected plane, say, 
x-y plane, and are computed with respect to the axes that 
are parallel and orthogonal to the line joining the two halos 
in the projected plane. Thus, the ellipticity component e\ 
(e.2) corresponds to the elongation and compression along 
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Figure 1. Ellipticity correlation functions of halos cn(r) as a 
function of the pair separation r. For the clarity, the results for 
different halo mass, from bottom to top, are multiplied by a fac- 
tor of 1, 10, 10 2 and 10 3 respectively, left panel — The results are 
presented for the halos of same spatial number density in the sim- 
ulations of 256 3 particles (open symbols connected with dotted 
lines) and in the simulations of 512 3 particles (solid symbols con- 
nected with solid symbols). From bottom to top, the halos have 
at least 20, 40, 80, and 160 particles respectively in the lower 
resolution simulations, and have 7 times more particles in the 
higher resolution ones (not 8 times because ag is slightly smaller 
in the latter), right panels — A few typical examples for the ellip- 
ticity correlation functions cn(r) measured in the simulations of 
512 3 particles in the boxes of 100 h — 1 Mpc (open triangles) and 
of SOOh^Mpc (solid circles). The halos more massive than Mh 
are included in the analysis, and the values of Mh (in Ii-'Mq) 
re indicated in the figure. The solid line is our fitting formula 
The resolution effect has been corrected according to the left 
panel. 



erally underestimated in the lower resolution simulations. 
The underestimation could be caused by an underestima- 
tion of the halo ellipticity and/or a poor determination of 
the halo orientations in the lower resolution simulations. We 
found that the latter might be the dominant cause, since the 
mean ellipticity of halos in the lower resolution simulations 
are found to be slightly higher (< 10%) than those in the 
higher resolution simulations. Nonetheless, the statistical re- 
sults converge very rapidly with the resolution of the halos, 
and iV m j n = 160 is sufficient for the simulation measure- 
ment. From this test, we found that the ellipticity correla- 
tion is underestimated by a factor of 2.0, 1.5, 1.25 and 1.05 
for iVmin = 20, 40, 80, and 160. These resolution effects are 
also confirmed when we compare the correlation functions 
from the simulations of different boxsizes with N p = 512 3 . 

A few examples for cn(r) measured in our simulations 
of N p = 512 are presented in the right panel of Fig. 1, 
illustrating for halos spanning over 2 orders of magnitude in 
mass. The shapes of cn(r) are very similar for different halo 
mass, and cn(r; Mh) for halos with a mass greater than 
Mh can be well fitted by 



en (r; > M h ) = 



3.6 x 10 _2 ( 



10 w h- 



'(YS^ + r 1 - 7 ) 



(2) 



which are shown in solid lines in the right panel of Fig.l, 
where r is in units of ft _1 Mpc. The results for the two dif- 
ferent boxsizes agree very well, as shown by the halos of 
M h = 6.4 x lO n /i _1 M0 (the second from top) in the figure. 
For the most massive halos M h = 5.2 x W 12 Ii~ 1 Mq, the cor- 
relation function appears to be slightly flatter than the fit- 
ting formula at small separation r < 1 /i _1 Mpc, though the 
discrepancy is small considering the errorbars of the simula- 
tion data. Here Mh is the mass of FOF groups with b = 0.1, 
which is about half of the nominal virial mass of b — 0.2 for 
typical CDM halos. 



(at 45° from) the line joining the two halos in the projected 
plane. For each sample of simulations, we compute the cor- 
relation functions for the three projections along the x~, y- 
and z-axes, and consider these projections as independent 
when estimating the errors. 

Our results confirmed that the cross-correlations C21 (r) 
and Ci2(r) vanish within the measurement errors (HRH00, 
CM00), and these functions will not be discussed anymore. 
As found in the previous studies (CM00), we find that the 
correlation £22(1") is very anisotropic, but cn(r) is nearly 
isotropic. For the presentation convenience, we will discuss 
the resolution effect in terms of cn(r). 

We analyze cn(r) for halos with particles more than 
iVmin in the simulations of 256 3 particles and of 512 3 par- 
ticles. The simulation box L is 100 /i _1 Mpc, and only the 
first three realizations are used for N p = 512 3 (so that the 
random phases of the initial fluctuations are the same in the 
two sets of the simulations). We consider halos with parti- 



cles more than JV n 



20, 40, 80, and 160 in the N p = 256 d 



simulations. Since erg is slightly lower in the N p = 512 3 sim- 
ulations, we fix the corresponding values of Amin for this 
set of simulations by requiring that the halos have the same 
number density n(> A m i n ), and we found AT m i n = 135, 270, 
540, and 1080. The results are compared in the left panel 
of Fig.l, which shows that the correlation function is gen- 



3 THE PROJECTED ELLIPTICITY 
CORRELATION FUNCTIONS 

The angular ellipticity correlation functions dj(9), as mea- 
sured in weak lensing surveys, are related to the three di- 
mensional correlations eij(r) through a modified Limber's 
equation (CM00, HRH00), 

c ((9) = J r 2 0(ri)ri^(r 2 )dridr2[l + ^12)}^ (r p , tt) 

J rj4>(ri)r%4>{r 2 )dridr 2 [l +f(ri 2 )] 

where n 2 is the comoving separation between the two galax- 
ies at n and r 2 , and r p and it are the comoving separations 
perpendicular to and along the line-of-sight. £(r) is the two- 
point correlation function, and <j>(r) is the selection function 
of a survey. For the flat Universe considered in this paper 
and at the small-angle limit (i.e. 9 <C 1 which holds well in 
large-scale lensing surveys), we have, 

ri + T2 n 

r p = — - — 6, ir = n — ri . 



(4) 



Since the correlations ejj(r) decrease rapidly at scales r ^ 
10/i _1 Mpc and the selection function <f>(r) changes much 
more gently on such scales, equation(^) can simplified to 

J rV 2 (r)drE lJ (r6») 



C«(0) 



[J r 2 4>(r)dr] 2 + J r 4 ^) 2 (r)dr J d7r£(r0,7r) 



(•») 
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Figure 2. The projected ellipticity correlation functions En(r p ) 
and T,22(r P ) (defined by eq.^j) of halos as a function of the project 
halo separation r p . The halos are selected and the results are 
presented in the same way as in the right panel of Fig. 1, but no 
vertical shifts are made for different halos in this plot. The solid 
curves are given by fitting formulae ("Q) and (^) respectively for 
the left and right panels. 
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Figure 3. Predictions for the angular ellipticity correlation func- 
tions for halos of mass larger than 1.2 xlO 10 , 9.6 xlO 10 , 6.4 X 10 11 , 
and 5.2 X 10 12 /i — x Mq (from bottom to top), assuming the source 
redshift distribution has a median redshift 1. For comparison, the 
observational data from deep lensing surveys of Van Waerbeke 
et al. (2000; open triangles) and of Wittman et al. (2000; solid 
symbols) are plotted. 



where we define Ejj(r 



dn[l + £(r p ,n)]cij(r p ,n) . 



(6) 



which we call the projected ellipticity correlation function. 
This function can be measured directly in the simulations. 
There are at least two advantages for adopting Eij(r p ). 
First, the three-dimensional correlations Cij(r) are generally 
anisotropic (in fact C22 is strongly anisotropic, see CM00), 
thus it is easier and more accurate to measure the projected 
functions. Second, more importantly, the projected elliptic- 
ity correlations contain all the information needed for pre- 
dicting the angular correlation functions dj (9) analytically 
(Eq. 0) once the redshift distribution of source galaxies is 
known. . 

We plot the projected ellipticity correlation functions 
£11 (r p ) and £22 (r p ) measured from our simulations of 512 3 
particles. The halos have correspondingly the same mass 
as those in the right panel of Fig. 1, and the resolution 
effects are corrected simply by multiplying the underesti- 
mation factors obtained from the convergence test (the left 
panel of Fig. 1). We find that the results from the two sim- 
ulation boxes agree very well. The functions En(r p ; ^ Mh) 
and E22(»>; ^ Mh) for halos with mass larger than Mh can 
be accurately fitted by the following expressions, 



En(r p ;> M h ) 



X 22 {r p ;>M h ) = 



r^ 5 (r p + 5) 
rg- B exp[i(*)3] 



■ ft -1 Mpc 



(7) 



■ ft~ 1 Mpc(8) 



where r p is in units of /i _1 Mpc. We note that from the figure, 
the readers may get an impression that our fitting formula 
underestimates E22 for the most massive halos (Mh = 5.2 x 
1O 12 /i _1 M0). But this impression is caused largely by the 



fact that the two simulation data points at r v — 1.2 and 
1.9 h~ Mpc, which are significantly smaller than the fitting 
formula, may be mistakingly considered as those of the less 
massive halos of M h = 6.4 x 10 n /i _1 M©. 

Van Waerbeke et al. (2000) and Wittman et al. (2000) 
have published their measurements for the angular correla- 
tion functions Cii(#) and C f 22(6 l ) from their deep surveys. 
Our projected functions can be used to predict the contri- 
butions of the intrinsic correlations to these surveys, if it is 
assumed that the galaxies have the same ellipticity as their 
host halos. We consider a redshift distribution function 



p(z) 



(-rexpR^f] 



(9) 



which approximately describes the deep surveys with a = 2, 
= 1.5 and z s w 0.7 (e.g. Smail et al. 1995; Wittman et 
al. 2000). Our predictions for the angular correlation func- 
tions, based on our fitting formulae for En(r p ;^ Mh) and 
E22(j>;JS Mh) and the fitting formula of Jing (1998) for 
£(r) of halos, are shown in Fig. 3, where halo mass Mh from 
lO lo /i _1 M to 5 x 10 12 /i _1 M Q is considered. If the source 
galaxies are more massive than 5X10 11 /i _1 M0, the contribu- 
tion of intrinsic ellipticity correlations Cn(ff) becomes com- 
parable to those observed. The intrinsic correlation C22{9) 
looks much smaller than the observation of Wittman et al. 
(2000), but the observed result is also significantly higher 
than the predictions of the LCDM model (cf. Wittman et 
al. 2000) for the lensing effect. We note that the intrinsic 
6*22(6') for M h > 5 x 10 11 h~ 1 M e is comparable to the pre- 
diction of the LCDM model for the weak lensing. 

In contrast to the previous studies of CM00 and HRH00, 
our results indicate that the intrinsic correlation of galaxy 
ellipticity can contribute significantly not only to shallow 
surveys but also to deep surveys. The difference of our con- 
clusions from theirs stems from two causes. First, CM00 
and HRH used the halos with more than 20 and 10 particles 
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respectively to predict the angular correlations, which can 
lead to an underestimation of at least a factor of 2 as our 
convergence test showed. Second, we show that the elliptic- 
ity correlation increases with the halo mass, and the galaxy 
mass can be higher than the mass 2.8 x lO'^^'Mg adopted 
by CMOO. 



4 CO NCLUSIONS 



We use p a set of high-resolution cosmological N-body sim- 
ulations to study the intrinsic ellipticity correlation func- 
tions of halos. With the simulations of different resolutions, 
we found that the ellipticity correlations converge once the 
halos have more than 160 members. For halos with fewer 



membe - s, the correlations are underestimated and the un- 



derestimation amounts to a factor of 2 when the halos have 
only 20 particles. After correcting for the resolution effects. 



we four d that the ellipticity correlations of the halos in the 
bigger box (L = 300 7l Mpc) agree very well with those 
obtained in the small box (L = 100 ft _1 Mpc). Combining 
these results from the different simulation boxes, we have 
presented accurate fitting formulae for the ellipticity cor- 
relation functions cn(r) and for the projected correlation 
functions En(r p ) and T,22(r p ) over a large range of halo 
mass (at least for 10 10 «S M h < 10 13 /i _1 M Q ). The latter two 
functions can be used to predict the contribution of the in- 
trinsic correlations to deep lensing surveys. With reasonable 
assumptions for the redshift distribution of galaxies and for 
the mass of galaxies, we found that the intrinsic ellipticity 
correlation can contribute significantly not only to shallow 
surveys but also to deep surveys, if the galaxies have the 
same shapes as their host halos. 
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